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In this paper we investigate the relationship between glassy and ferromagnetic phases in disordered Ising ferro- 
magnets in the presence of transverse magnetic fields, V. Iterative mean field simulations probe the free energy 
landscape and suggest the existence of a glass transition as a function of F which is distinct from the Curie 
temperature. New experimental field-cooled and zero-field-cooled data on LiHo a; Yi_ I F4 provide support for 
our theoretical picture. Here as well we present a collection of theoretical predictions for future experiments. 



The disordered Ising ferromagnet LiHoz Yi_2;F4 in a trans- 
verse field r provides a unique opportunity to study coexis- 
tent glassy and ferromagnetic tendencies. Application of T 
tunes both the glassiness and the ferromagnetism and, thereby, 
makes it possible to probe the complex free energy landscape. 
While the glassy characteristics of the paramagnetic trans- 
verse Ising systems have been studied both experimentally [[!]] 
and theoretically [B- [irj], as has the ferromagnetic transition 
in the pristine systern~[pl|, little is known about the interplay 
of ferromagnetism and glassiness. Equally interesting is an- 
other unique opportunity afforded by the LiHo a; Yi_ a ;F4 sys- 
tem, namely the ability to explore quantum mechanical effects 
arising from r ^ 0, in the presence of glassiness [|T2||. 

The goal of this paper is to study these two coexistent phe- 
nomena ( glassiness and long range magnetic order) via an ex- 
ploration of the free energy surface, and in this way, make sys- 
tematic predictions for the irreversibility characteristics (i.e., 
hysteresis, remanent magnetization, etc.). Ferromagnetism in 
the presence of significant disorder introduces a new class of 
challenging problems. We show here that the ferromagnetic 
onset temperature T c is distinct from the glass transition tem- 
perature T g , and that the latter is itself dependent on the his- 
tory of the system. Here, in contrast to the rather extensive 
theoretical studies (of the paramagnetic phase) in the litera- 
ture [S-|lQ], we associated T g , in more physical terms, with 
the onset of irreversibility in measurable characteristics. In 
general, there are three state variables, including longitudinal 
fields, H z , as well as T and T, which can be cycled in many 
different and non-commuting ways to arrive at a given mini- 
mum in the free energy surfaces. 

To support this physical picture, in this paper we also 
present new experimental data showing the difference be- 
tween field cooled (FC) and zero field cooled (ZFC) mag- 
netic susceptibilities as a function of T at fixed T in 
LiHoo.44Yo.56F4. We find that this difference vanishes at a 
(glass transition) temperature distinct from the Curie point for 
ferromagnetic ordering. 

We base our theoretical analysis on earlier work in conven- 
tional spin glasses [|l3]| and random field systems [|l4| which 
addressed the evolution of the free energy landscape using 
an iterative numerical mean field scheme, in which the re- 
action terms (which led to problems with numerical conver- 
gence [p3[]) were ignored. By following a given free energy 



minimum as it evolved with field and T, one arrived at rather 
good agreement between theory and experiment for the vari- 
ous history dependent magnetizations. 

This theoretical approach probes the system on intermedi- 
ate time scales, which are long compared to the time needed to 
"re-equilibrate" after a given free energy minimum has disap- 
peared (with temperature or field cycling), but short compared 
to the time needed to tunnel [fl|Jt2|] between meta-stable states 
on the free energy surface. 

The transverse Ising ferroglass is described by the Hamil- 
tonian: 

n = -\ E J v s ' s i - H * E s ? - r E s *> (1) 

where the sum (ij) is over the nearest neighbors, and the 
exchange coupling Jj j is given by the Gaussian distribution 
P{Jij) = v /l/27rJ 2 exp(-(J ij - J ) 2 /2J 2 ), where J is the 
variance and Jo is the shift. Here T oc Hf for small H t , where 
H t is the transverse magnetic field applied in the laboratory. 
We can obtain a mean-field equation for the average magneti- 
zation nii — (S*) T ({-) T denotes thermal averaging) at each 
site for S= 1/2: 

V 

mi = — tanh(/3£i), (2) 

where T t = E<i,j> Jij™>j + H z and ft = VW+^/ 2 - 
This mean field equation corresponds to minimizing the free 
energy F as a function of the set of mj. Here, the transverse 
field r effectively enters only through the modified Brillouin 
function of Eq. (^), reflecting the fact that T rotates the lo- 
cal spin axis away from the Ising or z direction. Introducing 
this term in effect mixes the "z-component" eigenstates of the 
r = problem. Finally, it should be noted that while the ex- 
periments we will present here address the magnetic suscepti- 
bility {x — J^Xi w ith Xi = drrii/dH z ), the calculations are 
based on the actual magnetizations. Iterative convergence of 
the susceptibilities at a given site has not yet been established. 

The system in the present study is composed of N x N x N 
spins with random bond configurations. Most of our exam- 
ples are for N = 20 and fixed Jo, although we studied larger 
(N = 40) size systems (and variable Jo) to verify conver- 
gence of our results. In order to make a closer connection to 
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FIG. 1. Magnetic field and temperature cycling protocols in the 
parameter space used to tune the free energy surface. The ZFC state 
is achieved through cooling from A to B and FC from A' to B'. 
The remanent magnetizations (called IRM and TRM) are obtained 
through the paths A-B-C-B, and D-C-B, respectively. The solid lines 
represent a consolidation of Fig. ^|as discussed in the text below. 

the experiments Jl2|], we present here the study of Jo = 0.2 J 
so that the system is only slightly ferromagnetic. For each 
set of parameters (r, H Z ,T), we start our iterations at the m% 
corresponding to the minimum of F evaluated at the previous 
T, H z , or r. We then update m, by solving the mean-field 
equations, Eq. (||), at each site until convergence is obtained at 
the nth iteration defined by 2~Zi( TO " ~ m i ) 2 /Ylii 171 ?) 2 — 
1CP 6 . Unless indicated, the magnetizations discussed in this 
paper are taken to be along the z- (Ising) axis. Finally, a small 
H z = 0.01J, was applied in all studies of glassy properties. 
This was needed to establish a fixed direction for spontaneous 
broken symmetry. 

Figure |l] is a three dimensional plot of the parameter space 
which we consider; these parameters cause the free energy 
surface to evolve in distinct ways. Various pathways in this 
parameter space will be referred to throughout the text. The 
solid lines shown here are irreversibility contours derived 
from magnetic hysteresis curves that are discussed later in the 
paper. 

In Fig. U we present an H z =0 phase diagram for the 
Hamiltonian of Eq.(l), omitting the lowest T regime where 
tunneling effects are important. There are two distinct lines 
separating different phases: the outer line indicates the phase 
boundary for the paramagnetic to ferromagnetic phases and 
the inner line is the ferromagnet to ferroglass phase boundary. 
The glassy boundary was determined by the condition that the 
transverse field cooled (FC) and zero field cooled (ZFC) mag- 
netizations are equivalent. The various cooling pathways can 
be seen in Fig. [j]. It should be stressed that in vanishing trans- 
verse field the quantity T g has meaning only as an asymptote, 
since the distinction between FC and ZFC becomes meaning- 
less when there is no field. The para-ferro boundary is defined 
by the set of (r, T) at which a spontaneous magnetization ap- 
pears. For a given N, this magnetization may be calibrated 




FIG. 2. Mean field phase diagram of the Ising ferroglass in trans- 
verse magnetic field T. Here PM represents paramagnet; FM, ferro- 
magnet; and G, glassy state. T g is determined by the point at which 
the FC and ZFC magnetizations are equal, and, consequently, is de- 
fined only for finite T. 

by first establishing a baseline zero, which is the magnitude 
of "spontaneous magnetization" due to finite size effects, es- 
timated from the Jo = case. By changing Jo, we were able 
to change the position of the para-ferromagnetic line relative 
to the glassy line; for smaller Jo, the para-ferromagnetic line 
will be inside the glassy line. 

We now turn to experiments. LiHoF4 is a three- 
dimensional, dipolar-coupled Ising magnet. In the classical 
limit ( r = 0), the Ho dipoles order ferromagnetically at 
a Curie temperature T c — 1.53K. Experiments confirm that 
standard mean field theory describes fully the critical be- 
havior of the phase boundary between paramagnet and or- 
dered ferromagnet jTllJ . Moreover, magnetically inert yttrium 
can be substituted for the holmium spins in single crystals 
of LiHo^Yi-zF^ permitting carefully controlled studies of 
the effects of quenched disorder. We suspended a needle of 
LiHoo.44Yo.5eF4 with T C (T = 0) = 670 mK from the mix- 
ing chamber of a helium dilution refrigerator into the bore of 
an 8 T superconducting magnet aligned perpendicular to the 
Ising c-axis (within 0.5°). A trim coil oriented along the Ising 
axis nulled stray longitudinal fields from the magnets. We 
measured the AC magnetic susceptibility \' m the frequency 
and excitation independent limits using a digital lock-in tech- 
nique. 

We compare in Fig. |3| the system response after various 
trajectories in H t — T space, where the static longitudinal 
field has been carefully set to zero. In the FC protocol, the 
sample was cooled in H t = IT from the paramagnet into 
the ordered state at T = 175 mK and then warmed; in the 
ZFC protocol, H t was only applied after cooling to 175 mK. 
The ferromagnetic transition is marked by a peak in field- 
cooled susceptibility at T = 295 mK. By contrast, the FC and 
ZFC susceptibilities bifurcate deeper in the ordered state, at 
T = 240 mK. As is consistent with the theoretical mean field 
phase diagram of Fig. we find clearly distinct signatures for 
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FIG. 3. Temperature evolution of the AC magnetic suscepti- 
bility x' °f LiHoo.44Yo.56F4 at Ht — IT for the FC protocol 
and the difference in %' under FC and ZFC protocols. The Curie 
temperature(T c ) in this transverse field is 295 mK while the FC and 
ZFC susceptibilities bifurcate at T = 240 mK. 

ferromagnetism and glassiness in the disordered Ising magnet 
LiHoo.44Yo.56F4. 

To elucidate these data, FC and ZFC processes were sim- 
ulated numerically as in classical spin glasses [OJ, with the 
fields applied in the transverse direction. We first cooled down 
the system in zero external field (path A — > B in Fig. [l]). 
We then applied a transverse field, and warmed up in tem- 
perature. In this way we obtained a zero-field-cooled magne- 
tization M ZFC as a function of temperature. These results 
are shown by the open symbols in Fig. [I|(a). We next cooled 
the system in the presence of a transverse field to T = 0.1J 
(path A' — > B') and then warmed up in the presence of this 
field. We thereby obtained the M FC curve (indicated by the 
solid symbols in Fig. f|(a)). The two magnetizations merge at a 
given temperature Tg, which we identify as the glass transition 
temperature of Fig. g. We found that M FC is reversible upon 
cooling and warming, while M ZFC is not reversible for a sub- 
sequent cooling from a temperature lower than T g . In this way 
we can regard the FC state as that which is closer to true ther- 
modynamical equilibrium, and it should not be surprising that 
this state has the larger magnetization of the two; this larger 
M takes advantage of the net ferromagnetic bias Jo. As the 
system is warmed the magnetization decreases monotonically. 
For both M FC and M ZFC , the magnetizations at the lowest 
temperatures decrease as V increases. This is a consequence 
of the off-diagonal components introduced by T, which act to 
reduce the net magnetization in the z direction. 

We turn now to the issue of remanent magnetizations. 
These remanences arise following the removal of a magnetic 
field, when the system becomes trapped in a metastable min- 
imum. In canonical Ising spin glasses there are two charac- 
teristic remanences [|l3]] (associated with H z ). These are the 
isothermal remanent magnetization (IRM) and the thermal re- 
manent magnetization (TRM). In the presence of a fixed trans- 
verse field we modify their definitions slightly, as follows. We 
first consider the case r = 0. When the system is cooled 
through path A — » B in Fig. |l] and a longitudinal field is ap- 
plied instantaneously (B — ► C) and subsequently adiabati- 




FIG. 4. (a) Temperature dependence of zero-field-cooled (ZFC, 
open symbols) and field-cooled (FC, closed symbols) magnetiza- 
tions for the Ising ferroglass at different T (b) Field dependence of 
IRM(open symbols) and TRM(closed symbols) versus H z for differ- 
ent T at T = 0.15 J. 

cally removed (C — > B), we refer to the resulting magneti- 
zation as the thermal remanent magnetization. When the sys- 
tem is cooled through path D — > C in a longitudinal field 
and the field is adiabatically turned off (C — » B), we refer to 
the resulting magnetization as the isothermal remanent mag- 
netization. For the present purposes we will be interested in 
the IRM and TRM at finite V. This corresponds to the paths 
A' -> B' -> C -► B' and D' -> B' -> C respectively. 
The behavior of the two remanences at different V is shown 
in Fig. @(b) for T = 0.15 J. Both the TRM (closed symbols) 
and the IRM (open symbols) become smaller as T increases, 
as expected since the spins are aligned more toward the trans- 
verse field direction. We find that the H z value where the two 
remanences are equal becomes smaller as T increases, indi- 
cating that at higher T it takes less H z to destroy the multiple 
minima on the free energy surface. 

As a final protocol, we address the behavior of magnetic 
hysteresis loops (for the magnetization as a function of vary- 
ing H z ), in the presence of fixed transverse fields. The results 
are shown in Figs. |5](a)-(d), where the transverse fields are 
r = 0.0 J, 1.0 J, 2.0 J and 3.0 J respectively. The loops were 
obtained by slowly decreasing the parallel field (H z ) from the 
high field limit (4 J) until an equally large negative field was 
reached and then sweeping back to close the loop. As is evi- 
dent, a non-zero transverse field enables the hysteresis loop to 
close at smaller H z . 

We may summarize these numerical hysteresis studies in a 
three dimensional plot of the associated irreversibility surface 
(solid lines in Fig. [I]). The surface consists of the locus of 
points in terms of the coordinates (T, H z , T) below which the 
system exhibits irreversibility, as reflected in magnetic hys- 
teresis loops. It should be stressed that this protocol is differ- 
ent from that used to obtain T g in Fig. [j]. While both protocols 
are experimentally accessible, the point at which irreversibil- 
ity sets in for a given (T,H Z ,T) is not unique, and, itself, 
depends upon the pathway to the point in question. This sce- 
nario is to be distinguished from the more conventional sit- 
uation which in which either T or H z is strictly set to zero, 
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FIG. 5. Hysteresis loops at fixed F =: (a)0.0J, (b)1.0J, (c)2.0J 
and (d)3.0J. Here T = 0.15 J. 



as in Refs [[jj]] and Refs [0] - JTo|], respectively. An impor- 
tant conclusion from Fig. l[is that raising the temperature or 
applying a longitudinal or transverse magnetic field progres- 
sively removes minima from the free energy surface until at 
sufficiently high temperatures or fields there is a unique state. 

The calculations which underlie this figure ignore thermal, 
as well as quantum, fluctuations. Given the latter, the low 
temperature regime should be viewed as inaccessible. Our re- 
sults can be contrasted with a recent low temperature study (in 
the paramagnetic phase) which investigated glassy phase in 
a quantum p-spin spherical model within an equilibrium sta- 
tistical mechanical approach Jlq], Here it was observed that 
hysteretic effects may also arise from a first-order phase tran- 
sition, rather than from the glassiness which we have empha- 
sized here. It should be noted that the experimental hysteresis 
reported in Fig. |3| is not strictly in the quantum regime. How- 
ever, to make a firm distinction between these two theoretical 
scenarios will require further experiments. 

In summary, in this paper we have emphasized the concept 
of history dependent measurements in transverse Ising ferro- 
glasses, a concept that has been widely recognized in other 
magnetic glasses Jl7l], Because of its focus on the complex 
free energy landscape, our approach should be contrasted with 
alternatives in the literature, which have also addressed the 
phase diagram of LiHo K Yi_ a; F4, but with an emphasis on the 
low temperature quantum regime using equilibrium statisti- 
cal mechanics. It should also be stressed that these transverse 
field configurations represent a unique opportunity to simulta- 
neously tune glassiness along with long range ferromagnetic 
order. Our predictions for the intermediate time scale behav- 
ior of the various history dependent magnetizations appear 
consistent with the FC and ZFC susceptibility measurements 
presented here, but further experiments will be needed to con- 
firm the sets of the predictions presented here. 
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